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ABSTRACT 



Photometric redshifts can be routinely obtained to accuracies of better than 0.1 
in Az/(l+z). However, the issue of dust extinction is one that has still not been well 

■ quantified. In this paper the success of two template-fitting photometric redshift codes 
00 \ (ImpZ and HYPERZ) at reliably returning Ay in addition to redshift is explored. 

. New data on the CNOC2 spectroscopic sample of 0.2<z<0.7 galaxies are presented. 

These data allow us to estimate Ay values from the observed Balmer decrements. We 
also investigate whether the empirical valu e of 7 = .44, th e ratio between gas- and 
star- derived extinction, as determined bv ICalzettil l)2001b|) . is necessarily the best 
value for this sample. 

When comparing the two codes to the Balmer-derived Ay, a correlation between 
Q-f the photometrically derived Ay, Phot- Ay, and the Balmer- Ay is found. The correla- 

q tion is improved when the empirical value of 7 = 0.44 is allowed to vary. From least- 

S_h ' squares-fitting the minimum in the reduced x 2 distribution is found for 7 ~ 0.25 ±0.2. 

^ \ For the sample of galaxies here, the factor of two difference in covering factor implied 

& ■ by the Calzetti ratio is found to be plausible. The CNOC2 galaxies with detected 

Balmer lines have some preference for an increased covering factor difference, which 
would perhaps imply they are undergoing more rapid, 'bursty' star formation than 

■ the galaxies Calzetti used in her derivation. 
5_i ' 

Key words: galaxies: evolution - galaxies:photometry - quasars: general - cosmology: 
observations 



1 INTRODUCTION 

Star formation rates (SFRs) and their global history (SFH) 
form the backbone of a slew of methods (observational, nu- 
merical, and analytical) investigating the processes of galaxy 
formation and evolution over cosmic time. The SFH is im- 
portant in indicating likely eras of dominant activity and 
in providing a self-consistent picture of chemical enrich- 
ment. This can then be compared to the predictions of 
semi-analytical models of galaxy formation and with inter- 
galactic medium (IGM) absorption line diagnostics. How- 
ever, SFRs can be imprecise due to complications arising 
from dust extinction. 
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COBE measurements of the cosmic far -IR/sub-mm 
background energy density iPueet et alJll996l) showed it to 
be e qual to, or greater than, the UV/optical background 
(e.g. lHauser et alJ Il998h implying that a large fraction 
of the energy from stars over the history of the Universe 
is hidden in the optical due to dust. The role of dust 
in high redshi ft galaxies has been discussed by many au- 
thors (e.g. Rowan-Robinson ct al. 119971: iPettini et aljll998l : 
ICalzetti fc HeckmarJll999l lAdelberger fc Steidell200Clh For 
example, star forming galaxies detected via the Lyman 
break technique at z~2 - 4 are estimated to be highly 
extincted in the rest-frame UV, meaning star formation 
rate s are ~ 3 — 10 times higher than if dust is ignored 
fe.g. lMeurer. Heckman fc Cal zetti 1999). Correction factors 
found for other high redshift star forming galaxies are of a 
similar order. The exact form of this extinction correction 
remains uncertain and in particular so does its evolution 
with epoch. 

Hence an important improvement for optical-based 



2 T.S.R Babbedge, R. Whitaker and S. Morris 



SFH studies is the determination of the extinction of galax- 
ies, in addition to their redshifts. In order to fully allow 
for variation from galaxy to galaxy, extinction needs to be 
measured as an a dditional free parameter to r e dshift . How- 
ever the study of iBolzonella. Miralles fc Pelld J2000l) found 
that the inclusion of Av as a free parameter in photometric 
redshift codes caused significant increases in alias i ng. I n 
a similar technique developed in Rowan- Robinson (2003), 
hereafter RR03, and extended in iBahbedge et a l. (2004), 
hereafter B04, these aliasing problems were reduced by 
setting several Ay priors. 

In this paper a sample of galaxies from the Cana- 
dian Network for Observ ational Cosmolog y (CNOC2) Field 
Galaxy Redshift Survey ijYee et, al]l2000f) is used to inves- 
tigate the reliability of Av values as determined by two 
SEP template fitting ph otometric redshift codes - HYPERZ 
(Bolz onella et al1l2000F ) and ImpZ (see B04) - by compar- 
ing the returned [z P hot, Ay] to the spectroscopically derived 
redshifts and Balmer decrement derived Ay's as calculated 
from CFHT MOS spectroscopic data. 

In fj2] the CNOC2 galaxy sample is set out, along with the 
follow-up CFHT spectroscopic data and Balmer decrement 
calculations. In i|21we discuss the link between Balmer ex- 
tinction and photometry. The photometric redshift method 
is briefly outlined in 2] and the results of applying the two 
redshift codes to the CNOC2 sample are presented in 33 
Overall discussions and conclusions are presented in JS] 

Note that for these investigations the flat, f2A=0.70 cos- 
mological model with Ho=72 km s _1 Mpc~ 1 is used. 



2 CNOC2 OVERVIEW 

The 2nd Canadian Network for Observational Cosmology 
Field Galaxy Redshift Survey (CNOC2) was conducted over 
a series of 53 night s on the Canada -France-Hawaii Telescope 
from 1995 - 1998 <Yee et alj EoOOV The survey covers a to- 
tal area of 1.5 sq. deg. spread over 4 patches equally spaced 
in RA. The dataset includes ~ 6000 galaxy spectra with 
spectroscopic redshifts to a nominal limit of R c ~ 21.5 in 
addition to 5-colour (J c , R c , V, B, U) 1 photometry of 40,000 
galaxies, complete to _R C =23.0 mag. The mean limiting mag- 
nitudes (Vega, 5cr) in each filter for the four areas are given 
in TableQand the filter response curves are shown in Figure 

in 

Use of a band-limiting filter restricted the spectral win- 
dow in the survey to 4400 - 6300A. This includes the [Oil] 
emission line for redshifts between 0.18 and 0.69. The sur- 
vey's resul ts on galaxy clustering over this range hav e been 
published jShenherd et al.l200ltlCa,rlher g et aljfennfft whilst 
an analysis of cosmic star formation history is in preparation 
(Whitaker et al, in prep.). 



2.1 Extension of the CNOC2 Survey 

In this paper we present follow-up observations of a subset 
of the CNOC2 sample of galaxies. The goal of the additional 



CN0C2 bands 




P 0.05 



2.0x10 3 4.0x10 3 6.0x10 3 8.0x10 3 1 .0x1 4 1 .2x1 4 
wavelength, Angstroms 

Figure 1. CNOC2 passbands combined with CCD response. 
From short to long wavelength they are U, B, V, Rc, Ic- 



Table 1. 5-sigma limits (Vega) for the four CNOC2 areas, and 
the Galactic extinction E(B — V) in each area. This can be con- 
verted to a correction in each band via Cardelli et al. (1989) with 
Ry=3.1. 



Area 


Ic 


Rc 


V 


B 


U 


E(B - V) 


0223 


22.97 


24.02 


23.95 


24.55 


22.98 


0.036 


1447 


23.52 


23.72 


24.35 


24.76 


23.27 


0.029 


0920 


22.85 


24.03 


23.94 


24.55 


23.18 


0.012 


2148 


22.69 


23.88 


23.70 


24.48 


23.08 


0.035 



observations was to observe wavelength ranges including the 
Ha and H/3 lines in order to set limits on the reddening of 
the objects and hence to allow estimation of the unreddened 
star formation rates for this subsample. As will be shown 
below, this subsample then allows us to test other means of 
reddening estimation and hence obtain reddening measure- 
ments for the entire CNOC2 sample. 

In order to maximise the run efficiency, masks were 
designed with slits assigned as a first priority to objects 
with 0.2<z spec <0.37 (allowing observation of both Balmer 
lines) and detected [Oil] A3727 emission. Second priority 
was assigned to objects within the same redshift range, but 
no detected [Oil], and third priority to galaxies within the 
CNOC2 magnitude range which did not yet have redshifts. 
This last sample will not be used in this paper. 

Data was taken w ith the CFHT MOS spectrograph 
llCrampton et al1ll992fl during a 4 night observing run in 
August 1999. 17 masks containing a total of 719 slits were 
observed, spread across 3 of the CNOC2 'patches'. The R300 
grism was used giving a potential wavelength coverage from 
4000-10, 000A (depending on slit location), with a dispersion 
of 5A per CCD pixel. A slit width of 1.5" was used giving 
a nominal spectral line FWHM of 3.4 pi xels or 17A. 

The data were reduced in IRAF jTodvl Il993l) 2 . The 



V photometry is actually g-band data calibrated to the V sys- 
tem based on Landolt standards 



2 IRAF, the Image Reduction and Analysis Facility, is a general 
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Table 2. Windows used for Ha and H/3 in this study 



Feature 


Window (A) 






continua 


line 


Ha 


6490-6536; 6594-6640 


6537-6593 


H/3 


4785-4815; 4911-4931 


4821-4901 



reduction was based on a slig htly modified ver sion of the 
CNOC2 standard reductions of lYee et ail feOOdh The steps 
included object finding, tracing and extraction, wavelength 
calibration, flux calibration, and interpolation over bad sky 
subtraction and zero order residuals. Error vectors were car- 
ried through the same procedure. 

2.2 Measurement of the Balmer Decrement from 
CNOC2 

Our own purpose writ ten code (W hitaker et al., in prep, 
based on the code by lBaloghl H999) was used to measure 
the fluxes and equivalent widths of the Balmer lines in each 
spectrum. In brief, the code measures the flux and equiv- 
alent width over pre-defined spectral windows, set for each 
line. Each window includes two continuum regions and a line 
region. We perform a 1.5a clip on the continuum regions to 
reject outliers and improve the quality of the continuum fit. 
Table [5] shows the windows we use for Ha and H/3 in this 
study. We no t e tha t our Ha definition is similar to that of 
iBaloeh et alJ J2002) with two minor modifications: Firstly 
we use a wider line region in order to fully encompass the 
Ha+[NII] flux in our measurement; this necessarily reduces 
the size of the blue continuum region by lA. This should 
have minimal affect on the quality of the continuum fits. The 
H/3 wi ndow we use is identical to that of lDressler fc Shecter I 
<1987ft . 

While it is true that the continuum fit is obtained from 
only 7 (4) independent points for Ha (H/3), the actual num- 
ber of data points used in the fit is considerably more: 23 
(13) for Ha (H/3), i.e. ~ 3 data points per resolution ele- 
ment. As a result, any spurious points, from cosmic rays or 
the like will be safely rejected without substantially reducing 
the quality of the fit. 

Figure [5] shows 3 typical spectra from our sample - the 
window regions are overplotted. 



Error estimates on the e quivalent widt h s are computed using 
the formulae set out bv lBohlin et alJ |[L983). For the error 
on the flux measurement we employ the following rule: 



where F and W signify the measured flux and equiva- 
lent width and of and aw their respective errors. 

purpose software system for the reduction and analysis of astro- 
nomical data. It is written and supported by programmers in the 
Data Products Program of the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



The original CNOC2 dataset contains multiple obser- 
vations of certain galaxies and using these we are able to as- 
sess whether we systematically under- or over-estimate the 
equivalent width errors. Our procedure is to compute for 
each pair of repeats a quantity e such that: 

e= m ~ W \ (2) 
(af + a? 2 )2 

Where W and a are as in Equation Q the subscripts 1 
and 2 refer to the separate observations. The perfect case is 
where the distribution of e is a unit gaussian of mean zero. 
This would imply that we have accurately measured the er- 
rors and there are no systematic under- or over-estimates. 

Such a unit gaussian is not found using the raw errors 
and so they are scaled according to the following equation. 

°true — Aa 2 aw (l + B(a 2 aw - ave 2 )) (3) 

Where A and B are multiplying factors to be deter- 
mined and ave is the average raw [Oil] equivalent width 
error larger than 2A. We recompute the distribution of e for 
a range of A and B and check whether it is a unit gaussian. 
The process is repeated until a unit gaussian is found. 

The same scaling we apply to the [Oil] errors is then 
applied t o the Ha erro rs in this work (see Whitaker et al., in 
prep., or lBaloghlll99fll for the full details of this procedure). 

2.3 Correction For [Nil] Emission and Stellar 
Absorption 

The window we employ in our spectral measurement code 
encompasses both the Ha and [Nil] lines. We cannot sep- 
arate the two lines in our data, nor are we able to accu- 
rately de-blend them using gaussian fitting tools (e.g. splot 
in IRAF) since the data is not of high enough resolution. 
Instead, t o correct for [Nil ] we assume an [NII]/Ha ratio of 
0.5 (as per lKennicuttll992ft . This is a simple approximation, 
however a by-eye examination of spectra in our sample re- 
veals that it is not obviously incorrect (see e.g. Figure 0. 
Using an extreme value for the NII/Ha ratio such as 0.33 
which lKennicutt fc Kent! jl983T) find for HII regions, results 
in a difference of ~ 0.4 in Av compared to the [NII]/Ha=0.5 
case. This value for HII regions does not incorporate the ef- 
fects of interstellar gas with its higher mean [Nil] /Ha ratio, 
and since our spectra are from the galaxy as a whole rather 
than individual HII regions we assume the 0.5 factor to be 
the most reasonable. 

The Balmer emission lines sit on top of stellar absorp- 
tion due to the presence of young and intermediate age stars 
in the line-emitting galaxy, so any measurement of their 
fluxes must take this into account. This is particularly true 
for H/3, where one often finds an emission line sitting in an 
absorption trough (see the middle right-hand panel of Figure 
[5] for an example of this). 

Again, due to the restrictions imposed by the data, we 
are not able to reliably fit the emission and absorption com- 
ponents separately within each line, so instead choose to 
apply the c orrections used by t he Sloan Digital Sky Sur- 
vey (SDSS: iHopkins et alJl2003T) . They measure a median 
absorption at Ha of 2.6A(note we do not employ their cor- 
rection of 1.3Asince we are using a window measurement 
method to obtain our fluxes, not the gaussian fits of the 
SDSS pipeline). For the H/3 correction, we use the value for 
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Figure 2. Three typical spectra from which we measure the Balmer decrement. Left column shows the Ha region, the right column 
shows the H/3. The window regions are overplotted: the solid (blue) lines show the continua regions, the dashed (red) lines the line region. 
The spectra are shown in rest-frame wavelength, the mean redshift of the whole sample is ~ 0.286 



Sb Galaxies (2 A) given bv lMiller fc Owenl |5p02), also used 
by the SDSS. 

In addition to the above two corrections to the flux mea- 
surements we also do the following: For sources with a secure 
Ha detection but either a low significance H/3 detection in 
emission, or detection at any significance level of H/3 in ab- 
sorption, we reset the H/3 flux to be its 3a error value. For 
the emission cases this is an obvious step, for the absorption 
cases however, it should be stated that the maximum H/3 
emission flux is still the 3a value - the absorption may of 
course be larger than this, but since the galaxies are emitting 
at Ha, they must therefore be emitting at H/3 also. 

Figures [3] and 2] show the flux distribution of Ha and 
H/3 both before and after the above corrections for [Nil] and 
stellar absorption. In figure^Jthe reset H/3 values are plotted 
as arrows, indicating an upper limit on the flux. 

2.4 Balmer Line Detections and Limits 

We divide the 719 objects observed in the extension to the 
CNOC2 dataset as follows: 

• 367 objects do not have redshifts in the original CNOC2 
survey and so are removed from our sample. 

• Of the remaining 352 objects, 46 are removed, either 
due to skylines obscuring the H/3 line (33 objects) , or prob- 
lems with the noise measurements arising from the data re- 
duction procedures (13 objects). 

This leaves 306 objects, and of this subset: 

• 46 objects possess both emission lines measurable at 
above 3a accuracy (we refer to these as 'measure' cases). 

• 153 objects have only Ha measurable to an accuracy 
above 3a - H/3 is not measurable to such accuracy in these 
spectra (we refer to these as 'limit' cases, i.e. we obtain a 
lower limit to the reddening for these objects). 



Flux Distribution Before Absorption and [Nil] Correction 
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Figure 3. The flux distribution of Ha and H/3 prior to corrections 
for [Nil] and stellar absorption. 



• 107 objects have neither Ha or H/3 emission detectable 
at above 3a accuracy, thus we are unable to estimate a lower 
limit on the reddening from these galaxies. 

Of the set of 46 measurements, five of them fall > 3a 
below zero and of these five, four are removed (see below). 
This leaves 42 measured sources for the analyses carried out 
later in this paper. We note that in a large enough sample, 
some objects are expected to lie 3a below zero reddening, 
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Flux Distribution After Absorption and [Nil] Correction 
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Figure 4. The flux distribution of Ha and H/3 after corrections 
for [Nil] and stellar absorption. Crosses are sources with a de- 
tection of at least 3a in both Ho and H/3, arrows denote sources 
with an Ha detection of at least 3a but where the H/3 flux has 
been reset to a 3a limiting value. Sources in which neither line 
could be measured at above 3a are not plotted. 



due to the nature of the statistics. We now describe the five 
objects in question: 

• Three objects appear to be blends where two (or more) 
galaxies are co-extracted from the slit; it is impossible to 
derive a reddening measurement from these spectra. 

• One object has an error in the noise vector. 

• One object has a high quality spectrum but a relatively 
small [Nil] flux. Our methodology of using a coarse cor- 
rection for the [Nil] flux serves to push extra objects into 
the 3er-below-zero reddening range (in addition to those we 
naturally expect from gaussian statistics) . This object is re- 
tained in the sample. 

Of the set of 153 lower limit cases we note that 70 pro- 
duce limits below zero reddening. Again, given our error 
approach, this is not surprising. 

The mean redshift of the sample of 'limits' and 'mea- 
sures' (Figure^J is ~ 0.262 while that of every object (Figure 
is ~ 0.286. 



2.5 Selection Effects 

The 'measures' sample we define and use in this work is 
subject to certain selection effects. 

The observations are clearly most sensitive to galaxies 
which have strong emission lines and low to moderate Ay. 
Any galaxies with weak emission lines and moderate Av will 
be lost from our sample since the lines will have been dimin- 
ished to the extent where they are undetectable. This is also 
true for galaxies with large extinction - at least H/3, if not 
both Balmer lines, will have been diminished by too great 
an amount to be detectable. The best that can be accom- 



plished in that case is to estimate a limit on the extinction, 
as we have done here. 

The problem of skyline contamination of the H/3 line is 
one which all projects attempting reddening meaurements 
in external galaxies will suffer from. However this will not 
affect our results unless there exists a strong correlation of 
extinction with redshift such that certain extinctions occur 
at certain redshifts where the H/3 line is obscured. Given 
the narrow redshift range of the survey and the strength of 
correlation that would be required this is not anticipated to 
be an issue. 



2.6 Calculation of extinction from the Balmer 
ratio 

The Ha and H/3 Balmer lines provide a probe of extinction 
for the optical region of the spectrum - their un-e xtincted in- 
tensi ty ratio is well known from atomic physics llOsterbrockl 
Il989h and is 2.86 in typical nebula conditions (T=10,000K, 
= 10 2 - 10 4 cm ). Their ratio thus leads to the Balmer 
optical depth, r, and can therefore provide a measure of the 
dust cont ent due to the intern al reddening along the line of 
sight fe.g lCalzetti et afl ll994'). down to an optical depth of 

T ~ 1. 

The observed ratio, R (, s , of Hq/H/3 can be shown to be 
related to the reddening to t he gas, E(B — V) , according 
to the following equation fe.g lCalzetti et aflll996f) : 



E(B-V) S 



2.5[log(R obs ) - log(Rj„t)] 
k(H/3) - k(Ha) 



(4) 



where Rint is the intrinsic Hq/H/3 ratio of 2.86 and 
k(Ha) and k(Hfi) are the extinction coefficien ts at the wave- 
lengt hs of Ha and Hj3 respectively (from ICalzetti et alJ 

The Calzetti curve JCalzetti et alJll993l) was used to 
compute Av's from the Balmer decrement, with the lowest 
measured Av being -2.25 (this is the source previously men- 
tioned that has low [Nil] flux) and the highest 5.15. These 
values are typically accurate to 30%; the mean of the ra- 
tio (Av/ctAv) is 0.27. Figure |5] shows histograms of the 42 
Balmer- measured Av values as well as the 153 limit cases. 

The choice of extinction curve is one of the sources of 
uncertainty - the e xtinction curves of local galaxies do differ 
(e.g see review bv ICalzettil l2001ah . The largest differences 
are due to the strength of the feature at ~ 2175A, thought 
to result from graphite grains, amorphous carbon or PAH. 
However, for the Balmer lines, this feature is not a factor: 
The Balmer lines are sampling the longer wavelength region 
where the extinction curves have a more constant slope, so 
that differences are reduced. The main difference between 
the curves here can be parameterised by the chosen value 
of the ratio Ry , which is the tota l-to-selective extinction 
ratio at V band iCalzetti et alll993T) . In these investigations 
Rv=4.05 is used. 



3 BALMER EXTINCTION VS. PHOTOMETRY 

Are we able to compare photometry-derived extinction 
values, hereafter Phot-Av, to Balmer- Av values in order 
to better understand the accuracy and/or reliability of 
photometry-derived extinction measurements? 
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Figure 5. The distribution of Balmer decrements (top panel), 
Balmer-derived Av values (middle panel) and associated error 
(bottom panel). The solid line shows results for the sources mea- 
sured at above 3(7 accuracy, the dashed line the limit cases (in 
this case the decrements and Balmer-derived Av values are lower- 
limits). 



The crucial point here is that the Balmer decrement 
probes the extinction Av of the ionised gas in the source - 
the extinction of molecular clouds/HII regions (star-forming 
regions), E(B — V) g . In the SED template-fitting method 
where photometry in a number of bands is compared to a 
library of templates with variable internal Av, the resulting 
best-fit returns a template type, redshift and Ay. In this 
case the extinction is due to the ISM acting on the stellar 
continua within the galaxy as a whole. Thus Phot-Av probes 
the dust obscuration of the stellar continuum, E(B — V)„. 
This is more complicated than just the total amount of dust 
between the observer and the source, as would be the case for 
a single star, as folded into this is information on scattering 
and the geometrical distribution of dust within the galaxy. 

Since the Balmer and photometry methods are prob- 
ing different regimes of the galaxy extinction, they might 
well be expected to give differing results. In a galaxy such 
as starburst where there is a significant amount of heavily- 
obscured star formation the two regimes will be linked - the 
galaxy's overall emission is dominated by the young star for- 



mation. At the opposite end of the scale an elliptical galaxy 
will have little dust or ongoing SFR. In this case the two 
measures are un-correlated (though both measures are likely 
to be low). For intermediate-type galaxies the relationship 
falls somewhere between these two regimes. 

The following relationship between the colour excess of 
the stell ar continuum an d of the nebular em ission lines is 
given in lCalzettil j2001bl) and ICalzettll dl997l) and is as fol- 
lows: 

E{B - V) s = (0.44 ± Qm)E{B - V) g (5) 

This follows on from ICalzetti et all (|l994h where 39 star- 
burst and blue compact galaxies were used to derive an 
effective extinction law. From this, it was found that the 
difference between the optical depth to the Balmer lines, 
and to the underlying stellar continuum at the lines was ap- 
proximately a factor two. Calzetti's relationship is a purely 
empirical one but it has often, in the past, been assumed 
to yield the most appropriate reddening corrections for the 
integrated l ight of extended star-forming r egions or galax- 
ies (e.g ICalzetti et aljEoQpT: IWestera et al]l2004h . Here we 
start with this ratio of 0.44 but then consider other possible 
values in order to see how the results are affected. 



4 PHOTO-Z METHODOLOGY 

The main outcome of applying a redshift code is the best- 
fitting redshift, extinction and template SED of each source. 
In the template-fitting procedure, the observed galaxy mag- 
nitudes are converted for each i th photometric band into an 
apparent flux, f° bs . The observed fluxes can then be com- 
pared to a library of template (T) fluxes, f* empl (z,T,Av), as 
calculated by convolving the template SEDs with the filter 
response functions. The reduced , Xred, ls computed for 
each point in the hyper-cube of redshift /template/ Av flux 
values and the best-fitting solution is selected. 

4.1 The templates 

Here, we use six galaxy templates, presented in B04; E, Sab, 
Sbc, Scd, Sdm and starburst galaxies. These were genera ted 
via spectrophotometric synthesis (see iBerta et alJl20o3 for 
more on this procedure) of several Simple Stellar Popula- 
tions (SSPs), each weighted by a different SFR and extin- 
guished by a different amount of dust, and were designed 
to reproduce in more detail the empirical low-resolution 
templates of RR03. These SSPs have been computed with 
a Salpeter Initial Mass Function (IMF) betw een 0.15 and 
120 solar masses, adopting the IPicklesI dl99ct) spectral at- 
las and extending its atmospheres o utside its original range 
of wavelengths with iKuruczl dl993T) models from lOOOA t o 
50,000A, as described in iBressan. Granato fc Silval ((,1998). 
Nebular emission is added by means of case B HII region 
models computed th rough the ionisation code CLOUDY of 
iFerland et alJ <ll998l) . The adopted metallicity is solar. In 
addition to the galaxy templates, two AGN templates can 
be considered. However, for the CNOC2 sample used here 
it was found that the AGN templates did not provide the 
best fit to any of the sources, (as expected based on their 
spectra) so in this paper only the six galaxy templates are 
considered - they are plotted in Figure |S| 
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One important consideration for the investigations here, 
and template-fitting codes in general, is that the extinction 
output by the photometric redshift codes refer to the Phot- 
Ay values of the best-fitting solutions. However, these values 
are those of the variable Av which has been added to the 
template in question. In addition to this there is also the 
issue of the inherent Ay of the templates. This is prob- 
lematic in that the templates originate from empirical tem- 
plates drawn from observations, along with alterations found 
to optimise redshift solutions in previous studies. Their re- 
generation via spectrophotometric synthesis where the over- 
all template is re-produced using a small number of simple 
stellar populations (SSPs) is designed to give some insight 
into the underlying physics of each template. Hence an over- 
all Ay can be defined, based on the Ay's of the contributing 
SSPs. 

In B04 these inherent Ay values were given (ranging 
from for the elliptical to 0.74 for the starburst). However 
it is possible to reproduce the same 3 overall template using 
different proportions of SSPs with differing Ay contribu- 
tions. For example, the starburst template can be generated 
from young SSPs to give the 'same' overall template, but 
with a total Ay of either 0.05 or 0.74. Similarly, the Scd can 
be generated with an Ay of 1.6 or 0.27. Clearly, then, we 
should not take the inherent Ay's of the SSP solutions at 
face value (we note that in the case of fitting SSPs to spec- 
troscopic data most of these degeneracies can be expected 
to be broken). Empirically, we would expect the inherent 
Ay of the late-type galaxies (Sab, Sbc, Scd, Sdm) to have 
an extinction to the gas of perhaps Ay ~ 0.5 — 1.0, and 
starbursts to be more dusty (Ay > 1.0, say). 

These SSP degeneracies would imply, then, that pho- 
tometric data is not able to constrain the extinction to an 
accuracy of less than perhaps ±1 in Ay. If this is the case, 
then one would not expect a correlation to be found between 
the Phot-Ay solutions and the Balmer-derived Ay even if 
the two methods were probing the extinction to the same 
regions of the galaxy. The results in ^will show that this is 
overly pessimistic, thus implying that SSP fitting considers 
more possible combinations than actually occur - there are 
additional constraints imposed in reality due to the nature 
of star formation, fueling and feedback that such fits do not 
incorporate. Since the details of inherent Ay are undeter- 
mined, inclusion in the analysis would introduce an addi- 
tional free parameter. For this relatively small sample there 
is insufficient data to properly constrain this so here it is 
merely noted that inherent Ay is an additional complica- 
tion and that it is expected to act to increase the Phot-Ay 
estimates. If the issue of inherent Ay is a large one then we 
will expect it to swamp any Phot-Ay/Balmer-Ay correla- 
tion - [J5]will show that, in practice, this is not the case. 

4.2 The ImpZ code 

This code was presented in B04. At high redshift, 

iMassarotti. Iovino fc Buzzonil i200ll) have shown that cor- 
rect treatment of internal (to the galaxy) dust reddening 
(the Interstellar Medium) and IGM attenuation (between 

3 here, 'same' means from the point of view of flux through a set 
of broad-band filters 



the galaxy and observer) are the main factors in photo- 
metric redshift success. The effect of internal dust red- 
dening for each galaxy is al terable via fitting for A y, us- 
ing the reddening curve of ISavaee fc Mathi a (^79 |), and 
the ef fects of the IGM are incorporated as in lMadau et alJ 
( 1996). Observed fluxes are compared to template fluxes for 
0.01s;logio(l+z)<0.90, equivalent to 0.02<z^6.94. 

Galactic extinction in the CNOC 2 regions is low ; using 
the extinction-wavelength relation in lCardelli et alJ Jl98STl 
we can calculate Galactic extinction in each of the CNOC2 
bandpasses and correct for it (see Table 0, though the re- 
sulting corrections are negligible. 

For the CNOC2 sample survey the following parameters 
were used, adapted from investigations presented in B04: 

Templates: The six galaxy templates (E, Sab, Sbc, Scd, 
Sdm and Starburst), with IGM treatment and Galactic ex- 
tinction corrections. 

Internal extinction: Ay limits of -0.3 to 3.0 in the Ay 
freedom. Negative Ay was allowed since the inherent Ay of 
the templates is non-zero. For the elliptical template Ay can 
take the value zero only, since ellipticals are not expected to 
have significant internal extinction. 

Magnitude limits: B04 found it necessary to apply 
redshift -dependent absolute magnitude limits to exclude un- 
physical solutions (such as super- luminous sources at high 
redshift). Here, the same limits are used: Absolute magni- 
tude limits of [-22.5-21ogio (1+z)]<M_b <— 13.5 for galaxies. 

Prior: ImpZ usually applies a prior expectation that the 
probability of a given value of Ay declines as |Ay| moves 
away from 0. This is introduced by minimising Xred + Q A V 
rather than % 2 (a=2 here) . This use of a prior can be viewed 
as a weak implementation of Bayesian methods and was 
reached at based on degeneracy studies in B04 and RR03. 
For this investigation it was of interest to turn off this prior. 
This is because the Balmer-Ay values of the sample range 
from -2.25 to 5.15, or for the starlight (multiplying by 0.44) 
-1 to 2.3. Hence applying an Ay prior would naturally tend 
to limit the Ay solutions to values in a similar range, regard- 
less of the input Ay limits, and so some correlation could 
be achieved simply because the sample all have Balmer-Ay 
values that equate to the lower end of the Ay parameter 
space. 

Tabled shows the parameters used in the final setup for 
ImpZ and HYPERZ. 



4.3 HYPERZ 

HYPERZ is described in full in iBolzonella et alJ fcOOd) . 

Here we present an overview of the code and the param- 
eters we used for each run. The HYPERZ code also uses 
SED fitting in its determination of photometric redshifts and 
for this work we chose to use the same templates as ImpZ. 
However, for the investigations presented here, we choose to 
apply HYPERZ for a 'best-case' scenario where both the 
redshift and template type is constrained in order to opti- 
mise the resulting accuracy of the extinction values. In order 
to create such a 'best-case' scenario, we choose to remove one 
degree of freedom by constraining the photometric redshift 
solution to within (z spec /100) of the known spectroscopic 
redshift. Additionally, we choose not to use the elliptical 
template in the analysis presented since none of the galax- 
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Figure 6. The 6 galaxy templates used. Dashed lines show the original RR03 templates (offset for clarity), solid lines shows the SSP 
generated versions, along with extension into the Far-UV (sub-lOOOA) as discussed in B04. 



ies with detected Balmer lines are ellipticals based on visual 
inspection of the spectra. 



We do not correct for Galactic extinction when run- 
ning HYPERZ. This correction is negligible in comparison 
to the photometric errors (see Tabled for extinction across 
the CNOC2 patches) and would not affect the results. Any 
non-detections in a given band were replaced with a flux 
of zero, the error being the limiting flux in the band (see 
Table [T] for 5cr limits). Ideally, one would want to provide 
the measured flux in an aperture placed at the location of 
the object even when that measurement is fainter than the 
nominal flux-limit of the survey and the resulting value has 
very poor S:N. However, the CNOC2 catalogue used here 
has had such measurements replaced by a 'non-detection' 
flag. Thus the best treatment in the template-fitting proce- 
dure is to restrict solutions to those that predict the flux in 
the non-detected band to be at or below the flux limit in 
that band. 



The HYPERZ code was run in its standard form, 
with the flat, nA=0.70 cosmological model with Ho=72km 
s _1 Mpc _1 used elsewhere in this paper and variable Av in 
the range <Av < 3. Table shows the parameters used in 
the final setup for ImpZ and HYPERZ. 



5 RESULTS 

5.1 ImpZ redshifts 

First, it is of interest to see how successful the photometric 
redshifts are when there is no Ay freedom - that is, Av=0 in 
the solutions. The results of this are plotted in the left panel 
of Figure |7| for 'measure' cases (black squares) and 'limit' 
cases (red crosses). It is immediately clear that not allowing 
Ay freedom has caused many of the ImpZ solutions to be in- 
correct - the code has been forced to substitute (incorrectly) 
additional redshift in place of the reddening action of dust. 
Indeed, there are a total of 11 of the 42 'measure' and 70 of 
the 153 'limit' sources (48% of the sample) whose photomet- 
ric redshifts lie outside the log(l+z apec )± 0.1 boundaries. 

If we now run ImpZ in the optimal manner, as discussed 
in H4.2I so that now Ay freedom is allowed, the solutions dra- 
matically improve. ImpZ was highly successful at deriving 
photometric redshifts in close agreement with the spectro- 
scopic values for 40 of the 42 'measure' and 147 of the 153 
'limit' sources (95% of the sample). Note we consider red- 
shift solutions in the range [0:7]). There were 8 (2 'measure' 
and 6 'limit') sources that obtained solutions at incorrect 
redshifts due to degeneracies in [z, template, Ay] space. We 
measure the accuracy of the photometric redshifts via the 
rms scatter <r z , calculated as follows: 
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Tattle 3. Final parameters for the two photometric redshift codes. 



ImpZ 



HYPERZ 



Templates 
Mj limits 
A„ limts, A„ step 
Reddening law 
Galactic extinction correction? 
strained to z. 
Cosmology 



E, Sab, Sbc, Scd, Sdm, Sb 
[-22.5-22o 9 io(l + 2)] <M S <-13.5 
-0.3^A„ $3.0, step 0.1 
Cardelli 1989 
Yes 

One run unconstrained, one run constrained 
n =l, O\=0.70, H =72km s^Mpc" 1 



Sab, Sbc, Scd, Sdm, Sb 
not constrained (but set by z spe c) 
0.0^A„ <3.0, step 0.27 
Calzetti 2000 
No 
Yes 

Qo=l, H A =0.70, H =72km s^Mpc" 1 



with N being the number of sources with both spectroscopic 
redshifts and photometric redshifts. For this data, the rms, 
(j z =0.32 when including the 8 outliers, or 0.07 when they 
are excluded. The comparison of spectroscopic and ImpZ 
redshifts is plotted in the middle panel of Figure |7| 

ImpZ is therefore successful at returning accurate red- 
shifts. However, can it also provide a measure of the extinc- 
tion compatible with that implied by the Balmer decrement? 
Since 2 of the 42 'measure' sources obtain an incorrect pho- 
tometric redshift this also means their extinction values are 
likely to be incorrect. In order to remove this (small) source 
of error in the following comparisons to the Balmer-derived 
extinction, we now constrain the redshift range explored by 
ImpZ to lie within 0.05 in log(l+z p hot) of z spec (plotted in 
the right panel of Figure |7J. Now, good solutions are found 
for all 42 'measure' cases and 153 'limit' cases, with <j z =0.06. 
It is the resulting ImpZ Av values from this setup that are 
considered from now on in the investigation. 



5.2 Range of ImpZ Av allowed values 

Although we wish to explore the accuracy of the extinction 
output from redshift codes by comparing to a sample with 
Balmer-derived measurements, we can also obtain an inter- 
nal estimate of how well constrained the [z,template,Av] so- 
lution is from the reduced x 2 distribution. For the solution 
with the minimum \ 2 , xLm the question can be asked: 
'What range of Av produces a fit at or near the correct red- 
shift (within 0.05 of log[l+z spec ]), with a reduced x 2 within 
y 2 ■ +1 ? ' 

The results of asking this question of each source is il- 
lustrated in Figure |S] It can be seen that for the majority 
of 'measure' cases (black), the Av is quite well constrained, 
in most cases to within 0.3 or so in Ay. Note that the lines 
in this plot can be discontinuous - for example, object 150 
in the plot has a best-fitting Av of 1.9 but has reasonable 
solutions in the Av range -0.3 to 0.4 and 0.8 to 1.1 which 
arise from fitting two other templates to the source. It is 
of interest to note that whereas 57% of 'measure' cases do 
not have discontinuous solutions, this is only true for 43% 
of the 'limit' (cyan) cases. Thus, determining the Av via 
photometry for these sources is more problematic, just as 
it is via the Balmer ratio method. Considering the 'mea- 
sure' cases in more detail, there are two objects where the 
Av is poorly constrained (objects 83 and 159 in the Figure, 
evidenced by their long black lines). Their best-fitting Av 
values are, respectively, 2.3 and 1.7, making them the most 
heavily extincted of the ImpZ 'measure' cases. Comparison 
to their Balmer-derived Av (3 ± 1 and —1.6 ± 0.5 respec- 



tively) would support the result for the first source but the 
Balmer decrement for the second source would imply neg- 
ative, or zero extinction, disfavouring the ImpZ best-fitting 
result or resulting in the interpretation that this source is 
problematic. 

Figure [§] illustrates the width of Av parameter space 
that lies within +1 of Xmin by plotting the distribution of 
this 'width' value (here the 'width' is simply defined as the 
maximum allowed Av minus the minimum allowed Av). For 
'measure' cases (black) the distribution drops quite steeply 
with width such that more than half the sources have a 
width of 0.4 or less. There is then a slight tail made up 
primarily of sources which had discontinuous solutions (such 
as one template with low Av and another template with 
higher Av) whilst the two sources with poorly constrained 
Av can be seen as a spike towards the maximum width of 
3.3 (i.e the full -0.3 to 3 Av range). The distribution for 
the 'limit' cases is more clearly bimodal, with a similar set 
of reasonably well-constrained sources with an Ay-width of 
0.4 or less, but a much larger set of sources with poorly 
constrained Av- Again, this is likely to be due to the nature 
of these 'limit' sources for whom the dust extinction is hard 
to determine (either via the Balmer lines or photometry). 

This analysis suggests that the photometric redshift so- 
lution has an inherently low extinction precision (at least 
with five band photometry) , such that Av is only precise to 
perhaps 0.3 for most sources, and is poorly constrained for a 
small subset. Rather than defining this internally estimated 
width value as the error in the Phot-Av value we choose in- 
stead to take the opposite approach for the comparison with 
the Balmer-derived Av- We will take the Phot-Av at face 
value and use the supposed relation to the Balmer-derived 
Av to provide an external estimate of the precision of the 
Phot-Av measurements. 



5.3 ImpZ - comparison to Balmer 

We can test how well the Calzetti ratio holds by comparing 
the Phot-Ay to 0.44*Balmer-A v . The distribution of A v 
residuals for ImpZ is plotted in Figure 1101 It can be seen 
that there is quite a spread to the distribution, though it 
is broadly centred on zero. Based on the findings in i|5.2l on 
the precision of the Phot-Av solutions some of this spread 
can be expected to arise from this low precision. Some of it 
can also be attributed to the accuracy of the Balmer-derived 
Av, which is typically accurate to around 30%. 

Figure 1111 plots (purple squares) the ImpZ Av values 
and residuals a s a function of Balmer- Av (multiplied by the 
ICalzett^OOlbl factor of 0.44). Note that the ImpZ A v values 
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Figure 7. Photometric redshifts from ImpZ. Comparison of spectroscopic and iMPZ-derived redshifts for 'measure' cases (black 
squares) and 'limit' cases (red crosses): Left panel shows results when Ay freedom is not considered; middle panel shows results for 
unconstrained redshift and Ay space; right panel shows results for redshift space constrained to be within 0.05 in log(l+z p ^ ot ) of z Bpec , 
but Ay unconstrained. Dot-dashed lines denote an accuracy of 0.1 in log(l+z), a typical photometric redshift accuracy. 
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Figure 8. ImpZ Ay allowed values. The range of Ay parameter space for each source that provides a solution with a reduced x 2 
within x^ in +l and that is at or near the correct redshift (within 0.05 of log[l+z spec ]). 'Measure' sources are shown as black lines, 'limit' 
sources as cyan lines. The best Ay solution value is indicated as a black ('measure' case) or cyan ('limit' case) square. 



have been taken at face-value and have not had an error 
assigned to them, since we wish to derive an error based on 
the comparison to the Balmer-Ay's. 

It can be seen that the residuals are smallest for Phot- 
Ay values of around 0.5 to 1 and that the residuals increase 
as we move away from this region (to either higher Ay 
or lower/negative Ay). Hence, the correlation to Balmer- 
derived Ay appears to be best for sources of intermediate 
extinction. It is also clear that none of the sources that were 
calculated as having negative Ay based on their Balmer lines 
obtain similar Phot-Ay values, lending weight to the sup- 
position that the Balmer method has fallen down for these 
sources due to limitations in the technique. The source with 
the largest Balmer-Ay, of 5.15T0.5, is also the source with 



the largest residual (of sources with non-negative Ay). Being 
so extincted, it is likely that this object is quite extreme, so 
disagreement between the star- and gas-derived extinction 
measures is to be expected. 

Calculating the following statistics: 

AA7 = V (Balmer[Av] - Phot[A v ]) /N (7) 

and 

a\ v = {Balmer[A Y } - Phot[A Y ]f /N (8) 
and outliers are defined by 

\AA V \ > 0.5 (9) 
gives AAy = —0.2, <ta v =0.83 and an outlier fraction, 
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Figure 10. A v residuals for HYPERZ (left panel) and ImpZ (right panel). The residual is (0.44*Balmcr[A v ]-Phot[A v ]). 
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Figure 11. Left: Av results for HYPERZ (cyan triangles) and ImpZ (purple squares). Plot is the Balmer[Av] (multiplied by the 
Calzetti (2001) factor of 0.44) versus the photometrically-derived Ay • Dot-dashed lines denote residuals of 0.3, 0.5 and 0.7 in Ay. 
Errors are not defined for the Phot-Av values. Right: Av residuals for HYPERZ (cyan triangles) and ImpZ (purple squares). Plot is 
Balmer-Av (multiplied by the Calzetti 2001 factor of 0.44) versus the residual (0.44*Balmer[A v ]-Phot[Av]). 



r), of 50 per cent. These results are better than one would 
infer from the large extinction degeneracies seen in the SSP 
fitting which were discussed in »j4. II Thus, a relationship 
between Balmer-derived and photometry-derived extinction 
measures is obtainable, though not a strong one. 



5.4 HYPERZ - comparison to Balmer 

Applying HYPERZ in a 'best-case' configuration (constrain- 
ing the HYPERZ redshift solutions to the spectroscopic 
values and excluding the elliptical template) gives similar 
statistics of AAv = —0.3, cta v =0.84 and an outlier frac- 
tion, r), of 50 per cent. The distribution of Av residuals for 
HYPERZ is plotted in Figure Q2| As with ImpZ, the distri- 
bution is broad but reasonably centred on zero. 



Figure lTTl plots (cyan triangles) the HYPERZ Av values 
and residuals a s a function of Balmer-Av (multiplied by the 
ICalzetti2001bl factor of 0.44). The residuals are again corre- 
lated with the Balmer-Av, being smallest for Phot-Av val- 
ues of around 0.5 to 1, and the negative Balmer-Av sources 
are again in poor agreement. 

Thus, HYPERZ and ImpZ portray a similar correlation 
to the Balmer-Av, though the agreement is noisy. 



5.5 Comparison of ImpZ and HYPERZ 

As well as comparing the two photometric redshift code's ex- 
tinction outputs to the Balmer-derived values, it is instruc- 
tive to compare them to one another to see if they tend to 
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Figure 9. Histogram of the width of ImpZ Av allowed 
values. Distribution of the width of Ay parameter space (defined 
by the minimum and maximum allowed Av) that lies within +1 of 
Xmin' 'Measure' sources are shown as a black line, 'limit' sources 
as a cyan dotted line. 




0.0 0.5 1.0 1.5 2.0 2.5 
MPZ A u 



Figure 12. Comparison between ImpZ and HYPERZ: Av 

results for ImpZ versus those from HYPERZ. Solid line is exact 
agreement, dot-dashed lines arc residuals of 0.3 in Ay. Note that 
for plotting purposes the values have been randomly altered by 
up to 0.02 in the x and y directions in order to separate points 
with the same/ very similar values. 



agree on a similar extinction value for a given source. A plot 
of ImpZ-Av versus HYPERZ- A v is given on Figure WI\ 

This shows that the two codes are in reasonable agree- 
ment about the extinction of a given source. Twenty- five 
of the forty- two sources (60%) agree within < 0.4 in Ay, 
thirty-five sources (83%) agree within < 0.5. The main dif- 
ference appears to be for five sources where ImpZ gives a 
high value of A v > 0.8 whilst HYPERZ tends to return 
a smaller Av estimate. Two of these sources (this includes 
object 159 mentioned in >15.21 have Balmer decrements that 
imply negative, or zero extinction, favouring the HYPERZ 
result or the interpretation that the sources are problematic. 
Two others have intermediate Balmer- Av consistent with ei- 
ther code's results and one has a larger Balmer- Av (this is 
object 83 mentioned in <I5.2I thus favouring the ImpZ result. 

Calculating similar statistics as when comparing to the 
Balmer- Av, comparison between the two codes' Av values 
gives AAy = —0.09, cta v =0.53 and an outlier fraction, rj, of 
17 per cent. This internal consistency check between the two 
codes gives increased confidence in the photometric redshift 
template-fitting method as a technique to obtain extinction. 



5.6 The Calzetti ratio 

We now consid er the empirical ratio of 0.44 introduced by 
ICalzettJ il997t) and set out in Eqn.|^| We consider a range 
of other ratios, 7, so that 

E(B - V), = lE{B - V) g (10) 

Varying this ratio from 0.05< 7 < 1, a x 2 minimisation 
analysis is carried out in order to find the ratio that 
gives the best correlation between the Phot-Av and the 
Balmer- Av, weighted by the errors in the Balmer- Av 
measurements. The results of this analysis are shown in 
Figure[ni For both ImpZ and HYPERZ, the lowest reduced 
X 2 is reached for 7 ~ 0.25. Based on the change in 7 that 



Table 4. Statistical results for the different comparisons between 
Balmer- Ay and Phot-Av ■ 



AAv o"A v outlier fraction, rj 

ImpZ, 7=0.44 dl2 083 50% 

HYPERZ, 7=0.44 -0.3 0.84 50% 

ImpZ, 7=0.25 dX4 072 38% 

HYPERZ, 7=0.25 -0.5 0.76 57% 



ImpZ/HYPERZ comparison -0.09 0.53 17% 



increases the reduced x 2 by one, the la error on this value 
is 7 = 0.25 ±0.25 for ImpZ and 7 = 0.25 ± 0.1 for HYPERZ. 

The resulting statistical measures, <7a v arm the outlier 
fraction, rj, are also plotted as a function of 7 for the ImpZ 
and HYPERZ results (Figure [1]}. The left panel shows 
how the rms in the residual, <ta v , varies with 7. A clear 
minimum is seen at 7 ~ 0.15 to 0.35 for ImpZ results, and 
at around 0.2 to 0.4 for HYPERZ. A similar minimum is 
seen in the range 7 ~ 0.2 to 0.35 for ImpZ results when the 
outlier fraction, rj, is plotted against 7 in the right-hand 
panel. For HYPERZ, the minimum is at around 7 ~ 0.3 to 
0.45. 

This analysis suggests that for this sample, the Calzetti 
ratio of 0.44 is a reasonable choice for 7 within the accuracy 
of the method, though a value of ~ 0.25 is preferred. A 
choice of 7=0.25 gives the following statistics: 

For ImpZ; AAv = —0.4, (ta v =0.72 and an outlier frac- 
tion, 7j, of 38 per cent. 

For HYPERZ; AAv = -0.5, cr Av =0.76, rj, of 57 per 
cent. 

As before, the distribution of Av residuals is plotted 
(Figure 1151 . It can be seen that the distribution is more 
peaked, though offset from zero. 
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Figure 14. Left panel: cr\ v for ImpZ (solid line with crosses) and HYPERZ (dot-dashed blue line, diamonds) as a function of 7, the 
chosen ratio between photometrically-derived and Balmer ratio-derived extinction measures. The Calzetti value of 7 = 0.44 is indicated 
as a long-dashed line. Right panel: Percentage outliers for ImpZ (solid line with crosses) and HYPERZ (dot-dashed blue line, diamonds) 
as a function of 7, the chosen ratio between photometrically-derived and Balmer ratio-derived extinction measures. The Calzetti value 
of 7 = 0.44 is indicated as a long-dashed line. 
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Figure 13. \ 2 analysis. The reduced \ 2 f° r ImpZ (solid line 
with crosses) and HYPERZ (dot-dashed line, diamonds) as a 
function of 7, the chosen ratio between photometrically-derived 
and Balmer ratio-derived extinction measures. The \ 2 values that 
are 1 above the minimum in the two distributions are indicated 
by horizontal lines (ImpZ, dotted; HYPERZ, dashed). 



Figure [Pol plots the Phot- Ay values and residuals as a 
function of Balmer- Ay (multiplied by 7=0.25). It can be 
seen that for lower Balmer- Ay, the Balmer- Ay tends to un- 
derestimate the extinction in comparison to the Phot- Ay 
value. If the negative Balmer-Ay are excluded, then the re- 
maining sources with positive Balmer-Ay do tend to follow 
the line denoting agreement, albeit with large scatter. 

In Figure tTTI the sources with only a lower limit on their 
Balmer-derived extinction (that is, a minimum 3<r Ha detec- 
tion but only a limit on the H/3 line detection) are plotted 
in comparison to the extinction as derived from the pho- 



tometric redshift codes. Here, no ratio 7 is applied to the 
Balmer Ay. Instead straight lines indicating different ratios 
are overplotted. Since these are lower limits, sources need to 
lie on, or to the left of a line to imply consistency with that 
chosen ratio. It can be seen that these lower limit sources 
are more consistent with lower values of 7. Of the 153 such 
sources, 146 (95%) are consistent with the 7=0.25 line when 
considering ImpZ solutions (purple squares), but only 130 
(85%) are consistent when considering the HYPERZ solu- 
tions (cyan triangles). 



6 DISCUSSIONS AND CONCLUSIONS 

A number of redshift codes routinely output an Ay value 
in addition to the best-fitting redshift, but little has been 
done to investigate the reliability and/or accuracy of such 
extinction measures. The main reason for this lies in the aim 
of such codes - they have been developed and optimised in 
order to derive redshifts. However, as the field of photometric 
redshift derivation matures, it is useful to consider some of 
the other parameters that redshift solutions output. 

In this paper we have asked the question 'Can a 
photometric redshift code reliably determine dust extinc- 
tion?'. The short answer would be: 'Not to a great accuracy'. 

Using a sample with extinctions derived from Balmer 
flux-ratios the Ay values produced by two photometric red- 
shift codes, ImpZ and HYPERZ, have been compared to the 
Balmer-Ay values. 

First, it was demonstrated that the inclusion of Ay was 
crucial in order to obtain photometric redshifts of high ac- 
curacy and reliability, such that 95% of the ImpZ results 
agreed with the spectroscopic redshifts to better than 0.1 
in log[l+z]). Without the inclusion of Ay freedom, there 
was a systematic and incorrect offset to higher photometric 
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Figure 15. Ay residuals for HYPERZ (left panel) and ImpZ (right panel) with 7=0.25. The residual is now (0.25*Balmer[Ay] 
Phot [Ay]). 
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Figure 16. Left: Ay results for HYPERZ (cyan triangles) and ImpZ (purple squares). Plot is the Balmer[Ay] (multiplied by 7 = 0.25) 
versus the photometrically-derived Ay. Dot— dashed lines denote residuals of 0.3, 0.5 and 0.7 in Ay. Errors are not defined for the Phot- 
Ay values. Right: Ay residuals for HYPERZ (cyan triangles) and ImpZ (purple squares). Plot is Balmer-Ay (multiplied by 7 = 0.25) 
versus the residual (0.25*Balmer[Ay]-Phot[Ay]). 



redshifts with many more incorrect redshift solutions. The 
existence of some negative Ay solutions may be indicative 
of the need for a bluer template in the template set, or for 
the inclusion of some additional free parameter in the fits. 
As the most important feature for template-fitted photo- 
metric redshifts is the location and identification of 'breaks' 
in the SED, the inclusion of Ay freedom can be seen, at 
first-order, as a modifier of the template SED's slope with- 
out having a strong effect on the breaks themselves. Hence 
a similar improvement may be achievable via a 'tilting' pa- 
rameter or similar which would act to alter the slope of the 
template SEDs. Since addition of dust extinction has a phys- 
ical basis, however, this is a preferable parameter, as long as 
we can demonstrate that there is some correlation between 



the best-fitting phot- Ay and the actual (or in this case, that 
measured via the Balmer ratio) dust extinction of the source. 
Thus, once the ability to derive good redshifts for the sample 
had been demonstrated a comparison between the phot- Ay 
and Balmer-Ay was carried out: 

The correlation between the Phot- Ay and the Balmer- 
Ay was similar for both codes, but was both noisy and not 
particularly strong. Based on direct comparison between the 
two codes, and investigations into the \ 2 solution space, a 
good part of this noise is derived from the inherent lack 
of precision the Phot-Ay solution has (perhaps 0.3 in Ay 
say), no doubt since it is based on only five photometric 
measurements. Additional noise arises from the precision of 
the Balmer-Ay, typically accurate to perhaps 30%, which 
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Figure 17. 'Limit' cases: Av results for HYPERZ (cyan triangles) and ImpZ (purple squares). Plot is the Balmer-derived Av 'limits' 
versus the photometrically-derived Av- Solid black lines denote 7 values of 0.25, 0.44, 0.6 and 1.0 for the conversion factor relating 
Balmer- Ay and Phot- Ay - Sources need to lie on, or to the left of a line to imply consistency with that chosen ratio. 



is due to the resolution of the spectrographs data. Given 
these errors, the correlation seen was in fact quite good. 

The correlation was improved somewhat when the 
empirical value of 7 = 0.44, the ratio between g as- and star- 
derived extinction, as determined bv ICalzettil (|2C)01b), was 
allowed to vary. From least-squares-fitting the minimum in 
the reduced x 2 distribution was found for 7 ~ 0.25 ± 0.2. 

The Calzetti ratio of 0.44 means that there is around a 
factor two difference in reddening such that the ionised gas 
(as measured by the Balmer decrement) is twice as reddened 
as t he stellar continua (as measured by th e photometry) 
(e.g iFanelli et alJll988t ICalzetti eral"lll994h . This implies 
that the covering factor of the dust is larger for the gas 
than for the stars, which can be explained by the fact that 
the ionising stars are short-lived and so for their lifetime 
remain relatively close to their (dusty) birthplace, whilst the 
majority of stars, contributing to the galaxy's overall optical 
luminosity are longer-lived and can migrate away from their 
dusty origins. 

For the sample of galaxies in this paper, this factor 
two difference in covering factor implied by the Calzetti 
ratio is found to be plausible, given the errors of the 
method. The sample has a some preference for an increased 
covering factor which would imply they are undergoing 
more rapid, 'bursty' star formation than the galaxies 
Calzetti used in her derivation. Perhaps more importantly, 
the results demonstrate the pitfalls of assuming that star- 
and gas-based extinction measures will give the same dust 
extinction given some conversion factor. Thus, correlation 
to Balmer-derived values are modulo the uncertainty in 
comparing star- and gas-based extinction measures. 



However the results presented here show that, given cer- 
tain considerations, there is potential in the application of 
photometric codes to reliably derive an extinction measure, 
though the precision is currently low. It is expected that the 
ability of photometric redshift codes to determine extinc- 
tion will improve with the availability of more photometric 
bands (here, there are five wide band filters between 3,000- 
9,000A). A sample with a combination of wide and narrow- 
band filters, with good wavelength coverage and range (in 
particular, extension to near-IR) will break many of the de- 
generacies and allow the codes to accurately differentiate 
between different possible fits. 

The results also show that it is important to note that 
this will be a measure of the star-based extinction, and will 
not necessarily be well correlated with the extinction to the 
ionised regions of a galaxy. 
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